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Abstract 

In this paper, we are studying the decoherence and entanglement properties for 
the two site Bose-Hubbard model in presence of a non-linear damping. We apply the 
techniques of thermo field dynamics and then use Hartree-Fock approximation to solve 
the corresponding master equation. The expectation values of the approximated field 
is computed self-consistently. We solve this master equation for a small time t so that 
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we get the analytical solution, there by we compute the decoherence and entanglement 
properties of the two-mode bosonic system. 

1 Introduction 

In recent years, there is lot of interest generated in the study of entanglement properties 
of ultra cold atoms [H El El IH El E]- In one such study, the single-site addressability in a 
two-dimensional optical lattice [S] has been demonstrated which could be a natural resource 
for applications of quantum information processing with neutral atoms. One of the popular 
model used to study the evolution of cold atoms and the BoseEinstein condensates in an 
optical lattice is Bose-Hubbard model [7]. In all the experimental demonstrations of ultra 
cold atoms, loss is an important role which gives rise to decoherence and in turn destroying 
the quantum correlations. In this paper, we examin the two-site BoseHubbard model to 
study the entanglement and decoherence properties of states under the action of non-linear 
damping. We consider the following master equation for density matrix p in a non-linear 
medium 



where a^, and are bosonic annihilation operators, k is a damping coefficient and H is the 
Hamiltonian for the Bose-Hubbard model which describes the optical lattice, discussed in 
the next section in detail . In ref O E] , the authors have studied the decoherence properties 
the Bose- Hubbard model with single mode damping. In this paper, we are studying the 
model in the presence of non-linear damping. 

For solving this master equation we use the techniques of thermo field dynamics and 
thereby the Hartree-Fock approximation to convert two-site Bose-Hubbard model into a 
two-mode bosonic system. The two-site Bose-Hubbard model is used to study Josephson 
tunneling between two Bose-Einstein condensates. The expectation values of the approxi- 
mated field is treated computed self-consistently. We solve this master equation for a small 
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time t so that we get the analytical solution, there by we compute the decoherence and 
entanglement properites of the two-mode bosonic system. 

The thermo field dynamics (TFD) [9l [TOl [HI [121 [13] is a finite temperature field theory. 
It is applied to many branches in high energy physics [10] and many-body systems [13] . The 
thermo field dynamics is used to solve the master equation by in presence of Kerr medium 
[in [15] using disentanglement theorem for any arbitrary initial conditions. This formalism 
presented in ref [lH [151 [El [13 [20] has two sailent features, first, solving of the master 
equation is reduced to solving a Schroedinger equation, thus all the techniques available 
to solve the Schroedinger equation are applicable here. Second, the thermal coherent state 
under the master equation evolution goes over to a thermal coherent state. 

A brief description of TFD is given in appendix D. In TFD, the master equation is given 

by 

= -zH\p) (2) 

where |p) is a vector in an extended Hilbert space ® Ti* and 

-iH = i{H ~H) + L. (3) 

where H is the Hamiltonian in Hilbert space H, H is the Hamiltonian in Hilbert space 
T-L* and L is the Liouville term, as mentioned in appendix D . Thus, —iH is tildian and the 
problem of solving master equation is reduced to solving a Schroedinger like equation namely 
eq ([2|). Then, the symmetry associated to the Hamiltonian such as su{l, 1) are exploited to 
solve master equation. 

In this paper, we follow this formalism to solve the master equation([TD for the two site 
Bose-Hubbard model in the presence of non-linear medium. 
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2 Two site Bose-Hubbard Model 



The Hamiltonian of the Bose-Hubbard model describes the interaction of bosons on a optical 
lattice and is given by 

k {k,l} 

+ V ^ + V ^ blbibkbk + ^Y. (^Wb\bi (4) 

^ k ^ k {k,i} 

with Ofc and bk bosonic annihilation operators referring to atoms in the internal states \Ni) 
and |iV2), respectively, with one boson in the kth lattice site and K is the number of lattice 
sites {A;, /} refers to the adjacent lattice points k and /. The interaction term J in the 
Hamiltonian describes the induced hopping between adjacent cells. The u is the frequency 
of the atom in the lattice. The on-site interactions of atoms are described by Ua and and a 
nearest-neighbour interaction by Uab- For further details see ref [7]. To study the decoherence 
and the entanglement properties of Bose-Hubbard model, for simplicity, we consider the toy 
model, in which the Bose-Hubbard model is written for the two site interaction only. The 
master equation ([1]) for the two site Bose-Hubbard Hamiltonian H (given in (jl])) is given by 
is 

d 

—p = —iuj{a) ap — pa) a) — iJ{a)bp — pa^b) — iJ{b'^ ap — pb'^ a) 
+i-^{a) a) aap — pa)a)aa) + i-^{b^b^bbp — pb^b^bb) 
—iuj{b^bp — pbb^) + i—^{a)b'^abp — pa)b'^ab) 

+ - {2abpa)b^ - a^b^abp - pa^b^abj . (5) 

At first we consider the special case to solve this master equation with J = Ua = Ub = 
and Uab = U, which corresponds to the Mott insulating phase. Then the master equation 
reduces to 

d U 

—p = —iuj{a'ap — pa^a) — iuj{b^bp — pbb^)+i — {a^b^abp — pa'^b^ab) 

+ — (2abpa^b^ — a^b^abp — pa^b^ab^ . (6) 
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Now, we apply the thermo field dynamics techniques to convert the master equation into 
a Schroedinger equation by applying |/) from the right to the eq 

||P) = -^^|P), (7) 
where |p) extended Hilbert space H "H* and the Hamiltonian is given by 

— iH = —iu{a^a — aa^)—iu{b''b — bb^)+i—(a^b^ab — aJ}a^b^) 

+ ^ (2abab - a'^b^ab - aba'^b'^) , (8) 

here the a, b are the annihilation operators act on the H and a and b are the annihilation 
operators act on the Hilbert space "H* ( for detail see appendix D). 

This master equation ([8]) is a non-linear equation, and in general it is difficult to get 
an analytical solution for this equation. A way out would be to apply the Hartree-Fock 
approximation and treat the approximated field self-consistent ly. In the case of thermo field 
dynamics, a selfconsistent theory using Hartree-Fock approximation is developed for the non 
linear master equations in presence of the non-linear medium in ref [20] ■ Applying, Hartree 
Fock approximation for each term in ([8]) as follows 

d)h^ab = {ab)^ abab = ab{ab), dbab = ab{ab), (9) 

with (ab) = (ab) = A(t), the Hamiltonian in (|H]) is decoupled into tildien and non tildian 

parts 



iH = i{H^ + H2), (10) 



where 



and 



2 
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The solution of is then given by 

|p(t)) = {exp[-i I dtHi] ® exp[-i I dtH2])\p{0)), (13) 



where |p(0)) is an initial state in H (S) H*. 

It is clear from the above that the two Hamiltonians Hi and H2 are independent in the 
sense that Hi is acting on non-tildian system and H2 is acting on tildian system Hence, we 
can work with one of the Hamiltonian and similar thing work for the other Hamiltonian 
(exept for interchanging between u and —u). 

To study the decoherence and entanglement properties of the two-mode states under the 
action of the master equation ([7]), we would like to exploit the underlying symmetry associ- 
ated with the Hamiltonians f lTT]) and f[T^ . This is accomplished by defining the following 

Af = a^a + b^b, /C+ = a^b\ /C_ = ab (14) 
= dd^ + bb\ IC+ = d^b\ t_ = db (15) 

which satisfy the 5^(1, 1) algebra 

[A/',/C+] = /C+, [A/',/C_] = /C_, [/C+,/C_] = 2Ar. (16) 

Similar algebra holds for tildians. Rewriting the Hamiltonian f lTT]) in terms of su{l, 1) gen- 
erators one gets 

-tHi = -^oo^^+^^{IC.-IC^)+^^^{IC^ + IC.), (17) 

and similarly the Hamiltonian ( !T2|) . It is clear that the Hamiltonians are associated with the 
sm(1, 1) symmetry. Hence, the underlying symmetry of the master equation ([7]) is sm(1, 1) x 
s-u(l,l). One of the important features of this sm(1, 1) symmetry is that it gives rise to 
squeezing and in turn it gives rise to entanglement. If the initial state is a two- mode Gaussian 
state then, under the action of two-mode squeezing, the final state turns out to be also a 
Gaussian state. Thus, one can use the separability criterion for Gaussian states [23] to check 
for the entanglement of the time evolved state. 
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3 Solution of Master Equation 

To get the solution to the master equation ([7]) we have used the Hartree-Fock approximation 
and treated the approximated field as background field, and the later has to be computed 
self-consistently in order to get the solution of (^^. For doing a self-consistent analysis we will 
first exactly diagonalize the Hamiltonians eqs ( 1TT|) and ( IT2|) with the help of the underlying 
su{l, 1) symmetry. It is evident form the eqs (11 II) and (|T2|) that the Hamiltonian in eq ([7]) is 
decoupled into tildien and non tildien parts as mentioned in eq (llOp . Hence, we work with 
one of the Hamiltonians and similar analysis goes through for the other Hamiltonian. By 
considering the Hamiltonian 

H, = a;(ata + 6t6) + !^M(^a5_at6t)_^^(at6t + ^5)^ (18) 
and applying the following transformation 

A = na + iy*b\ A^ = fi*a^ + ub (19) 
B = i2b + v*a\ B^ = ij,*b^ + ua (20) 

where fi = e*^|/i| and u = e*^|z/|, we diagonalize the Hamiltonian Hi. Similar analysis goes 
through for the other Hamiltonian H2. With a bit of algebra for n = 0, one can exactly 
diagonalize the Hamiltonian f|T8|) by using the transformations f|T9|) and fl20|) . Let k = 
for t = 0, then we evolve the master equation at this time. Then the corresponding non 
tildien Hamiltonian Ho is given by 

Ho = u{a^a + b^b)-^^^^{a^b^ + ab), (21) 

and the final diagonalized Hamiltonian Hf (which is diagonalized version of Hq), after the 
unitary transformation fll9p and (120]) . is given by 

Hj = S-^ {r)HoS{r) = [A^A + B^B + 1] (22) 

where 

S{r) = exp[r}C- — r*/C+] = exp[ra'^b'' — r*ab] (23) 
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here = Au"^ — f/^A^(0) and r is related to and v in eq (fT^ and (I2U|) via the following 
Bogolyubov coefficients : 

w , / X f/A(0) , , 

u = coshir) = I . „ , = sinh(r) = — , ^ ==. (24) 

Note that the above satisfy 

|;i|2 _ |iy|2 = 1. (25) 

This fixes A'^{0) = ^^j^ = Constant. The solution to the Schrodinger 

th^lMt)) = -iH^mt)) (26) 
where Hamiltonian Hq is given by eq fl2Tl) . is given by 

l^o(t)) = exp[r/C_ - r*/C+]|^/'o(0)) = exp[ra%^ - r*a6]|7/'o(0)). (27) 
Similar solution exists for the Hq 

Ho = uj{a^a + b^b)-^:^^^{d^b^ + ab). 



(28) 



3.1 Self Consistency Analysis 



Now, we compute the background field A(t) self-consistently by taking the initial state of the 
two-mode system to be the vacuum state. In thermo field dynamic notation the two-mode 
vacuum state is |p(0)) = |00,00) and in the usual notation p{t) = e~*/"'*-^|0, 0)(0, 0|e*/"'*-^. 
In thermo field dynamics the background field A{t) is computed as expectation value of the 
averaged creation and annihilation operators, and is given by 

A(t) = {ab) = {I\ab\p{t)) = Tr{abp{t)). (29) 

After doing a bit of algebra (see Appendix B) one gets A(t) to be 

A(t) = (1 + zcjt) A(0) + t dt'A{t')sinh\UA{Q)t). (30) 

2 JO 
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This is nothing but the Fredholm integral equation of the second kind. Computing the 
integral upto first order one gets 

A{t) = (1 + zut)A{0) - ^ '-[- - —sznhi2Qt)]. (31) 

By considering Qt to be small (so that sinh{2Qt) ~ 2Qt) one gets 

A(t) = (1 +iwt)A(0). (32) 

Hence the solution of master equation ([7j) is given by eq (fT3|) . which when written in terms 
of su{l, 1) generators becomes 

\pit)) = expHH'tMO)) 

= {exp[Ca3^f + Ca-IC- + Ca+JC+] ® explCbs^f + Cb-JC- + Cb+}C+]) |p(0)), (33) 

here Cas = i(^t, Ca- = I dt^{iU Ca+ = I dt^{iU-K), (bs = i^t, Cb- = I dt^{iU +h) 
and = / dt^Y'{iU — k). Using disentanglement formula [22] one can write eq fl33|) as 

\p{t)) = {{exp[Ta+IC+]exp[ln{Ta3J\f]exp[Ta-}C_]) 

® [exp[Tb+)C+]exp[ln(Tb3Af)]exp[n-}C-]) )|p(0)). (34) 

where 

_ 2Ci±sinh(j)i _ 1 

2(t)iCosh(t)i - Ci3Stnh(j)i (cosh(i), - ^^sinhcj),) 

with 

0? = f -G+O- (36) 

and i stand for a and h. 

It can be clearly seen that the Fj are functions of the background field A(t). By using 
the expression for A(t) from eq flH21) (which is valid for small t) we get 

r., . + (37) 



By taking [ill ± k) = —(e'^^ then one gets 



+ 



(38) 



2 ' A ' 

By considering any arbitrary initial state p(0) = Z^m « Pm,"!"^) (^1 of ^ single mode system 
where |m) and |n) are number of states, in the thermo field dynamic notation, there by using 
eq (17^ (of appendix D), |p(0)) takes the form 



(39) 



Putting eq ( 139)) in one gets 



min(m' ,n') oo 

E E 

g'=0 p'=0 

min{rn,n) oo 

X E E 

q=0 p=0 



V 

Y 
A 



m' + p' — q' 
pi 

m + p — q 
p 



^ n' +p' -q' 



V 



P' 

n + p — q 
p 



' m' \ 









n 



jp' ^p^^j{m'+n'-2g'+l)/2j 
X Pm+p+p'-(g+g'),n+P+p'-(g+9') (0) • 



(m+n-2ij+l)/2 |--p^ jg 



(40) 



4 Entanglement 

The solution of the master equation ([6]) is a pure state in thermo field dynamic notation and 
is given by 

|p(t)) = {exp[Ta+K:+]exp[ln{Ta3^f]exp[Ta^K:^] 

xexp[Th+}C+]exp[ln{T,3Af)]exp[T,^K:^])\p{0)). (41) 

where are given in eq f l38|) . In this master equation, as there is no mixing between the 
tildian and non-tildian modes in the Liouville space, one can write the solution in the system 
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Hilbert space as 



p{t) = (exp[r,+/C+]exp[/n(r,3Anea;p[r,_/C_]|^(0))(^(0))| 
{exp\ra+lC+]exp[ln{Ta^,N]exp\ra-K,J\) 



(42) 



where again Ti± are given in eq (155]) . One can clearly see that this a two-mode pure squeezed 
state of the two-mode Hilbert space. It is well known that two-mode squeezing gives rises 
to entanglement. As an example, we take the initial state to be two-mode thermal state. 
Then, to calculate the entanglement of the time evolved state we go over to phase space by 
following transformation 



1 

71 



1 

71 



X - ip^ 



(43) 
(44) 



Putting them in the eq f l42|) one gets the two-mode squeezed thermal state in the phase space 
corresponding to the two-mode the real squeezing transformation 



/ 



S{r) 



cosh{r) 


sinh{r) 






cosh{r) 


—sinh{r) 



sinh{r) 


cosh{r) 




\ 

—sinh{r) 


cosh{r) I 



(45) 



where the squeezing parameter is given by r = ^(1 + ^Xt. The covariance matrix of a 



two-mode thermal state is given by 



(46) 



where the ni and n2 are sympletic eigenvalues. The covariance matrix for the state p(t) in 
eq fH2|) is given by 

V = S{r)aS\r) (47) 
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V 



\ 



\ 



(4J 



where 

^ p -s 

p s 

-s g 

s g y 

with p = niCOsh'^{r) + n2sinh'^{r), q = nisinh'^{r) + n2C0sh'^{r) and s = -[- "-i+^z sinh{2r). 
So the covariance matrix is of the form 

A C 
B 



V 



(49) 



where 



A 













p 1 





q 
q 



C 



-s 
s 



J 




Then the separabihty condition [23J for any two-mode state reads as 
detAdetB - 



1 \ 2 1 

- - \detC\j - tr{AJCJBJC^J) > -{detA + detB) 



(50) 



(51) 



So, according to condition (15T|) in eq (H2|) . the state p(t) is entangled if and only if the 
condition is satisfied. Then the thermal state is entangled 



sinh (r) > 



(52) 



(ni + n2)2 

For rii = n2 = n, the amount of entanglement in p{t) is given in terms of logarithmic 
negativity 

E^{r) = -^[Log{e-'^/n)]. (53) 



5 Decoherence 



As we have seen the previous section, that the solution to p(t) of the master equation ([6]) in 
the Hilbert space H is given by 

pit) = iexp[TaMexp[lniTa3Af]exp[Ta-}C^])\^lj{0)){^{0)\ 

{exp[Ta+}C+]exp[ln{Ta3Af]exp[Ta-IC^]) (54) 
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where again Ti± are given in eq fl55]) . To calculate decoherence effects of p{t) we compute 
and is given as follows note that 



2+2 



Pit) = exp[-A(0)2(l + ^)Ct]p(0) 



(55) 



Then 



2+2 



tU t 

Tr[p'^{t)] =Tr[Y,{m,n\p'^{t)\m,n)] = exp[-A(0)2(l + ^)Ct] 



(56) 



The behaviour of decoherence is plotted fig. 2. One can see immediately that for the short 
time it self as the value of damping coefficient increases the system decoheres faster. 




Figure 1: Time vs Entanglement: Here u = .25, A(0) = .25 and n = .5 the blue, pink, 
orange and green corresponds to ( = 0.1, 0.25, 0.5 and 0.75. 



6 Conclusion 

In this paper, the techniques of thermo field dynamics and the Hartree-Fock approximation 
are used to solve the master equation for the special case of two-mode Bose-Hubbard model 
in the presence of a non-linear damping. We have treated the approximated field self- 
consistently and computed it for a small time t analytically. Then, the decoherence and the 
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entanglement for this system are computed. We show that the entanglement for this system 
for a short time increases when the initial state is in two-mode thermal state. We interpret 
this behaviour due to the existence of the non-linear medium. To get the exact picture for 
the behaviour of the entanglement for a long time, one has to do the numerical studies. It can 
be seen form the decoherecnce plot that as the value of the damping coefficient increases the 
damping in the system is faster, as expected. We expect that the further numerical studies 
using this model will give better results and these results can be applied to condensed matter 
systems. 



7 Appendix 



7.1 Appendix A: Calculation of Bogolyubov Coefficients 

Consider the following Hamiltonian 

Ho = u{a^a + b^b)-^^^^^{a^b^ + ab), 
By applying the following Bogolyubov or squeezing transformation 

A = na + u*b\ = ^*a) + vb 



(57) 



B = fib + u*a\ = fi*b^ + 



va 



(where = e^^\^\ and v = e*^|z/|), we diagonalize the Hamiltonian Hq. Let us first rewrite 
the Hamiltonian in a convenient form as 



Ho = W b] 



1 m 


n \ 


a 


\ n 


m j 


6t 



(5J 



where m = iu and n 
(1581) with 



-iUA{0) 



. By identifying the coefficient matrix appearing in the eq 



^ m n ^ 








^ n m ^ 





cosh{r) sinh{r) 
sinh{r) cosh{r) 



(59) 
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and with a suitable normalization \J\m\'^ — \n\'^, the Bogolyubov coefficients can be read off 
as follows 

f/A(0) 



jjL = cosh{r) 
The above satisfy 



^,2 C/2A2(0/ 



V = sinh{r) 



2^ 



2 _ ;72A2(0) 
4 



I |2 I |2 1 



(60) 



(61) 



7.2 Appendix B: Computation of A(t) 
Consider the following field variable appeared in eq fill I) and 0121) 

A(t) = (ab) 



(62) 



The initial state being the vacuum state the thermo field initial state is given by |p(0) 
|0,0,00) and in usual notation 

p(t) = e-*/"'*^|0,0)(0,0|e^/°'*^ where H is is equal to in eq (dH). So we have 

{ab) = {I\ab\pit)) =Tr{abp{t)) =Tr{abe-'I'^'^\0,0){0,0\e'I'^'^) 
= ^(m, n|a6e-^/°'*^|0, 0)(0, 0|e^/°'*^|m, n) 

m,n 

= 5](0,0|e^/'^*^|m,n)(m,n|a6e-*/'^*^|0,0) 

m,n 

= (0,0|e*/'^*^(^ \m,n){m,n\)abe~'I ^^'^10,0) 

m,n 

= (0,0|e'/'^*^a6e-'/'^*^|0,0) 

= {0,0\{l + i J dtH)ab{l-i J dtH)\0,0) 

= {0,0\{ab + i J dt[H,ab])\0,0) 

= {0,0\ab) +i J dt{0,0\[H,ab]\0,0) 

= PP + I y'rft(0,0|[Ca3A^ + Ca-/C_+Ca+/C+,/C_]|0,0) 

= ixu + i J dtuniy + J dt(a3\lJ'\'^ 
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= (1 + tiot)A{0) + J dtA{t)stnh\r) 

= (1 + tujt)A{0) + J dtA{t)stnh\UA{0)t). (63) 

7.3 Appendix C: Calculation of and F 

Considering only for small time t, i.e., neglecting A^(t) and higher order terms one has from 

m 



-^4 4 ' 

which gives 

20±g^n/i(^) ^ 20±3zn(f; _^ 

icut(cos/i(^) + stnh{'^)) ut ' ' ^ ' 

Here we have used cosh{x) + sinh{x) = and —isinh{ix) = sin{x). Again considering ut 
to be small, that is using sin{ujt) = ut one has from eqf l65|) : 

= =C.±(i-^) (66) 



Thus 



it^(,U±.)il-f). (67) 



Then putting the value of A(t) from eq eq fl^ into eq fl^71) one gets 

. %)(.f;±.)(, + 4),l_!^)l^„f;±.)(l + ^). (69) 
7.4 Appendix D: Thermo Field Dynamics 

A brief description of TFD is given below. The dissipative term in any master equations 
makes it difficult to apply the usual Schroedinger equation techniques with pure states to 
mixed states. The thermo field dynamics (TFD) provides such a formalism. In TFD the 
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mixed state averages are expressed as scalar products and the dynamics is given in terms of 
Schroedinger like equation. A density operator p = |A^)(A^| corresponding to a Fock state 
|A^) in the Hilbert space "H is viewed in TFD as a vector p = \N,N) in an extended Hilbert 
space H^H*. The central idea in TFD is to construct a density operator 1/2 < a < 1 
as a vector in the extended Hilbert space H<Si'H*. 

Here the averages of operators with respect to p reduces to a scalar product: 

{A)^Tr[Ap] = (/-"I^Ip"), (70) 

where is given by 

|p°) = p"|/),with, p^^p'^^I, (71) 
where |/) is the resolution of the identity 

in terms of a complete orthonormal basis {|n)}^g in H. The state vector |/) takes a 
normalized vector to an another normalized vector in the extended Hilbert space ® "H*. 
The matrix A{a, a^) acts like A /. 

(It may be noted that for any density operator the states |p"),l/2 < a < 1 have a 
finite norm in the extended Hilbert space 1-1®%*. This is not in general true for the state 
1/2 < a < 1, which includes |/). These state are regarded as formal but extremely 
useful constructs.) 

The creation and the annihilation operators a\a) , a, and a are introduced as follows 

a\n,rh) = \/n|n — l,m), o)\n,rh) = \/n + l|n + 1, m), (73) 
a\n,fh) — \/rn\n, m — 1), a}\n,fh) — \/rn + l\n,rn -]- 1) . (74) 

The operators a and o) commute with a and . It is clear from the above that a acts on 
the vector space % and a acts on vector space %* . From the expression for |/) in terms of 
the number states 

n 
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it follows that 

a\I) = aV), a^I) = a|/), (76) 

and hence for any operator A (written in terms of a a"^ and their complex conjugates), one 
has 

A\I)=A^\I), (77) 

where A is obtained from A by making the replacements tilde conjugation rules a — )■ a, — >■ 
, «—)■«*. An immediate consequence of this is that the state |p") which remains unchanged 
under the replacements a ^ d, ^ d'^ , and c number — )• complex conjugates C by applying 
the the hermiticity property of p i.e. = p. The tildian property reflects the hermiticity 
property of the density operator. 

The evolution of a conservative system in terms of is given by the von Neumann 
equation 

|p"W = ^[^,p1, (78) 

by applying |/) from the right and one gets 

d_ 
di 

where 



p"(t)) = -z^Ip"), (79) 



-iH = i{H - H). (80) 

In TFD, one can derive a Schroedinger like equation for any state |p") with arbitrary value 
of a. For dissipative systems the evolution equation is given by master equation 

-p(t) = ^{Hp-pH) + Lp, (81) 

where L is Liouville term. The non equilibrium thermo filed dynamics is developed and 
analysed in a = 1 representation. Hence, from now on, we work in, a = 1 representation 
fl9[ [20| [18]. In this representation, for any hermitian operator A, one has 

{A) = {I\A\p) = {A\p)=Tr{Ap). (82) 

18 



By applying |/) to the eq f ET]) from the right one goes over to TFD and the master equation 
is given by ([2]). with —iH being a tildian and thus problem of solving master equation is 
reduced to solving a Schroedinger like equation namely eq ([2]). 

Historically the thermo field dynamics was developed in a = | representation. In this 

1 

representation |po ) is related to the |0, 0) by a unitary transformation, which is nothing but 
the Caves-Schumaker state, for details ref P [IHl [IB [El [I3l [li EI] . 
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2 4 6 8 10 

Figure 2: Time vs Decoherence: Here u = .25 A(0) = .25 the blue, pink, orange and green 
corresponds to C = 1, 0.75, 0.5 and 0.25. 
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